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Phys.  Plasmas  3,  2098  (1996)]  is  used  to  accelerate  planar  plastic  foils  to  velocities  that  for  the 
first  time  reach  1000  km/s.  Collision  of  the  highly  accelerated  deuterated  polystyrene  foil  with  a 
stationary  target  produces  ~Gbar  shock  pressures  and  results  in  heating  of  the  foil  to  thermonuclear 
temperatures.  The  impact  conditions  are  diagnosed  using  DD  fusion  neutron  yield,  with  ~  106 
neutrons  produced  during  the  collision.  Time-of-flight  neutron  detectors  are  used  to  measure  the 
ion  temperature  upon  impact,  which  reaches  2  —  3  keV. 
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I.  INTRODUCTION 

Here  we  discuss  experiments  that  explore  the  acceler¬ 
ation  of  near  solid  density  planar  targets  to  extremely 
high  speeds  (~  1000  km/sec)  using  a  krypton  fluoride 
laser.  These  velocities  exceed  those  needed  for  even 
the  highest  velocity  designs  contemplated  for  inertial  fu¬ 
sion  implosions  (400-500  km/s)  [1].  Collisions  of  the 
high  velocity  targets  provide  a  means  to  obtain  pressures 
above  1  Gbar,  a  regime  of  interest  to  high  energy  den¬ 
sity  physics.  The  temperatures  and  pressures  are  high 
enough  at  1000  km/s  to  produce  thermonuclear  neutrons 
when  employing  targets  containing  deuterium.  The  ve¬ 
locities  reported  here  are  to  our  knowledge  the  highest 
achieved  by  laser  acceleration  of  material,  and  approach 
those  thought  to  be  needed  for  impact  ignition  [2]. 

Aside  from  using  a  specially  designed,  long  (tens  of 
km)  linear  accelerator  [3],  laser  driven  ablation  appears 
to  be  the  only  method  of  accelerating  bulk  matter  to 
such  speeds.  Speeds  of  700  km/s  have  been  achieved 
previously  on  a  third  harmonic  Nd:glass  (A  =  351  nm) 
Gekko/HIPER  laser  at  Osaka  University  [4].  There  are 
several  challenges  to  achieving  speeds  ~  1000  km/s.  It 
requires  coupling  ~  0.5  GJ/g  of  kinetic  energy  to  the 
target  while  keeping  the  target  near  solid  density,  or  at 
the  very  least,  above  the  critical  density.  The  speeds  re¬ 
quired  are  close  to  or  higher  than  the  exhaust  velocity, 
requiring  most  of  the  target  mass  to  be  ablated.  Further¬ 
more,  given  a  limited  laser  energy  the  target  is  necessar¬ 
ily  thin,  making  it  more  susceptible  to  heating  due  to 
radiation  transport  from  the  corona  and  to  breakup  due 
to  the  Raleigh-Taylor  instability  (RTI).  RTI  effectively 
limits  the  distance  over  which  the  target  can  be  acceler¬ 
ated  before  breaking  apart.  Achieving  the  highest  previ¬ 
ous  velocities  on  a  glass  laser  required  the  use  of  special 
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techniques  to  mitigate  RTI,  such  as  double  ablation  using 
high-Z  dopants  [4]  and  two-color  irradiation  [5].  A  way 
to  mitigate  this  situation  is  to  keep  the  non-uniformities 
seeding  the  instability  to  a  minimum  and  apply  increased 
pressure  to  accelerate  the  target  over  a  shorter  distance, 
while  making  use  of  ablative  stabilization  to  reduce  the 
RTI  growth  rates  [6] . 

Utilizing  a  short  wavelength  (248  nm)  KrF  laser  makes 
the  task  significantly  easier.  Shorter  wavelength  results 
in  ablation  at  higher  densities,  giving  higher  ablation 
pressure  for  a  given  intensity.  It  also  allows  higher  in¬ 
tensity  on  target  before  the  onset  of  deleterious  laser- 
plasma  instabilities  [7].  Finally,  the  application  of  in¬ 
duced  spatial  incoherence  (ISI)  smoothing  [8]  gives  an 
ultra- uniform  illumination  of  the  target,  limiting  the 
seeding  of  the  RTI. 

This  paper  is  organized  as  follows.  Sec.  II  describes 
the  experimental  setup,  Sec.  Ill  gives  the  experimental 
results,  Sec.  IV  discusses  the  results  together  with  an¬ 
alytic  modeling  and  simulations;  finally,  conclusions  are 
drawn  in  Sec.  V. 


II.  EXPERIMENTAL  SETUP 

Fig.  1  shows  the  experimental  setup.  Targets  consist 
of  two  foils:  the  first  one  10-16  /xm  thick  polystyrene, 
the  second  15-20  /xm  thick  graphite.  The  laser [8]  with 
the  pulseshape  and  spatial  profile  shown  in  Fig.  2  irradi¬ 
ates  the  first  foil,  causing  it  to  ablatively  accelerate  and 
impact  the  second  foil.  The  distance  between  foils  is  var¬ 
ied  between  500  and  700  /xm.  Both  foils  are  significantly 
wider  than  the  laser  spot  size  («  2  mm  on  the  side  vs. 
0.46  mm  FWHM,  respectively).  The  trajectory  of  the 
flyer  is  obtained  by  imaging  it  side-on  onto  the  slit  of 
a  streak  camera  using  monochromatic  x-ray  imaging  at 
1.86  keV  with  a  spherically  bent  quartz  crystal[9].  The 
image  is  either  from  backlighting  and  self-emission,  or 
self-emission  only  (no  backlighter).  The  amount  of  self- 
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FIG.  1:  (Color  online)  Schematic  of  the  experimental  setup. 
Shown  is  the  laser-illuminated  polystyrene  foil  followed  by  the 
graphite  witness  foil.  Also  shown  shematically  is  the  side-on 
imaging  using  Bragg  reflection  from  a  curved  crystal.  The 
adjustable  contrast  aperture  placed  in  the  image  plane  of  the 
backlighter  source  controls  the  amount  of  self-emission  x-rays 
reaching  the  detector.  The  lower  inset  on  the  right  shows 
an  image  of  a  shot  with  backlighting  and  self-emission.  The 
upper  inset  is  a  shot  with  self-emission  only. 
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FIG.  3:  Sample  side-on  streak  camera  image  with  backlighting 
and  self-emission.  The  laser  is  shining  from  the  left.  Visible 
are  shadows  of  the  two  foils  with  the  accelerated  part  of  the 
first  foil  in  between.  Strong  x-ray  self-emission  is  visible  at 
the  point  of  impact. 


III.  EXPERIMENTAL  RESULTS 

A.  Acceleration 


FIG.  2:  Laser  pulseshape  (left)  and  spatial  profile  lineout 
(right),  t  —  0  corresponds  to  the  half-rise  of  the  laser  intensity 
on  target. 


emission  imaged  is  controlled  by  an  aperture  placed  at 
the  image  plane  of  the  backlighter  source [9].  Backlight¬ 
ing  is  done  using  a  Si  target  irradiated  by  dedicated  laser 
beams.  Timing  of  the  streak  camera  with  respect  to  the 
laser  pulse  is  set  using  the  time  of  half-rise  of  x-ray  self¬ 
emission  from  the  target  front  surface,  with  an  estimated 
uncertainty  of  ±100  ps. 

The  impact  of  the  accelerated  first  foil  generates  a 
strong  shock  in  both  foils.  Visible  light  emission  from 
the  shock  breakout  on  the  rear  of  the  second  foil  is  im¬ 
aged  onto  a  slit  of  a  streak  camera  using  an  f/3  optic. 
The  polystyrene  is  pure  CH,  pure  CD,  or  either  of  the 
two  with  0.4%  (atomic)  Br  doping.  Five  neutron  time- 
of- flight  (nTOF)  detectors  are  positioned  at  various  an¬ 
gles  around  the  target:  three  at  approx.  2.6  m  from  the 
target,  and  two  at  3.6  m. 


A  sample  side-on  streak  camera  image  is  shown  in 
Fig.  3.  X-rays  from  both  backlighting  and  self-emission 
are  imaged  in  this  shot.  The  vertical  shadows  are  of 
the  first  foil  (12/xm  deuterated  polystyrene  doped  with 
0.4  at%  of  Br)  and  of  the  graphite  second  foil.  The 
shadow  of  the  accelerating  part  of  the  first  foil  is  visible 
in  between.  Self-emission  is  visible  initially  from  the  first 
foil  during  the  initial  peaked  part  of  the  pulse.  Strong 
self-emission  appears  at  the  point  of  impact  with  the  sec¬ 
ond  foil  starting  at  about  1.6  ns. 

Without  the  Br  doping,  the  trajectory  of  the  acceler¬ 
ated  target  is  not  so  readily  observable  with  backlight¬ 
ing,  in  particular  because  carbon  opacity  at  1.86  keV 
drops  dramatically  as  temperature  rises  above  100  eV. 
Self-emission  is  a  much  higher  contrast  technique  in  this 
case.  A  sample  image  of  self-emission  only  with  the  tra¬ 
jectory  traced  out  is  shown  in  Fig.  4.  In  this  case  the 
backlighting  is  turned  off  and  the  contrast  aperture  is 
fully  open  to  allow  all  the  self-emission  at  1.86  keV  re¬ 
flected  from  the  crystal  to  be  imaged  onto  the  streak 
camera.  In  this  way  the  self-emission  from  the  target  is 
visible  throughout  the  acceleration,  and  peaks  strongly 
when  the  impact  occurs. 

The  target  trajectory  can  readily  be  traced  in  the  im- 
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FIG.  4:  Sample  side-on  image  of  the  accelerating  10.5  /mi  CH 
target  taken  with  self-emission  only.  The  trajectory  traced 
from  the  image  is  overlaid  as  a  dashed  line.  The  impact  on 
graphite  marked  by  bright  emission  is  clearly  visible. 


FIG.  5:  Velocity  of  the  10.5  /mi  CH  target  obtained  from 
the  side-on  image  in  Fig  4.  The  acceleration  increases  toward 
the  end  as  the  target  mass  is  ablated  away.  The  error  bars 
(one  standard  deviation)  obtained  from  error  analysis  of  the 
trajectory  fit  are  shown  in  gray. 


age  as  the  boundary  of  self-emission  (shown  as  dashed 
line  overlaid  on  the  image  in  Fig.  4).  The  trajectory  is 
fitted  with  a  polynomial  and  differentiated  to  obtain  the 
target  velocity.  The  velocity  from  the  trajectory  of  Fig.  4 
is  shown  in  Fig.  5. 

The  procedure  used  to  fit  the  trajectory  and  obtain  the 
velocity  is  as  follows.  At  each  spatial  pixel  of  the  image, 
the  time  pixel  of  the  half-rise  of  the  self-emission  is  found. 
The  resulting  [£,  x{t)\  data  points  are  least-squares  fitted 
with  a  fourth-degree  polynomial.  Fourth  degree  polyno¬ 
mial  was  chosen  for  physics  reasons  as  well  as  because 
it  reduces  the  x2  error  of  the  fit  [10]  as  compared  to  a 
lower  degree  polynomial.  A  second  degree  polynomial 
would  imply  constant  acceleration  of  the  foil,  whereas  an 
increase  in  acceleration  is  expected  as  the  mass  is  ab¬ 
lated.  A  third  degree  polynomial  would  give  constant  in¬ 
crease  in  acceleration,  however  variation  in  acceleration 
is  expected  to  be  non-linear  due  to  the  decrease  in  laser 
intensity  in  time  as  well  as  target  mass  reduction  owing 
to  ablation.  Increasing  the  polynomial  degree  further  re¬ 
sults  in  higher  uncertainty  in  the  fit  without  a  physical 
reason  to  do  so. 

The  fourth  degree  polynomial  obtained  from  the  fit  of 
the  trajetory  is  then  differentiated  to  obtain  the  velocity 
as  a  third  degree  polynomial.  The  error  bars  on  velocity 
have  been  obtained  by  two  independent  methods.  The 
first  method  assumes  that  the  errors  in  the  trajectory 
points  are  normally  distributed  and  uses  the  matrix  of 
covariances  of  the  errors  of  fit  coefficients  of  the  trajec¬ 
tory  and  error  propagation  [11].  The  second  method  uses 


a  Monte  Carlo-like  “bootstrap”  method  [10],  where  the 
error  estimate  is  obtained  by  redoing  the  trajectory  fit 
many  times  on  data  sets  synthesized  from  the  original 
data  by  replacing  a  random  fraction  of  the  points  with 
others  from  the  same  data  set.  The  two  error  estimates 
agreed  to  better  than  20%.  The  larger  of  the  two  was 
used  for  the  one  standard  deviation  error  bars  shown  in 
Fig  5. 

The  velocity  trace  shows  the  acceleration  decreasing  as 
the  drive  intensity  drops  (Fig.  2),  then  increasing  again 
after  t  =  0.9  ns  as  the  drive  intensity  levels  off  and  as 
a  significant  portion  of  the  target  mass  is  ablated  away. 
The  increasing  acceleration  toward  the  end  is  also  promi¬ 
nent  in  the  simulations,  as  will  be  seen  in  Section  IV  B. 
The  target  is  moving  at  just  over  1000  km/s  prior  to  the 
impact. 

Fig.  6  shows  a  streak  image  of  the  visible  light  emit¬ 
ted  from  the  rear  surface  of  the  second  foil  for  the  same 
shot  as  Figs.  4  and  5.  Shock  breakout  appears  at  t  — 
1.24  =b  0.05  ns.  This  is  just  after  the  collision  seen  on  the 
side-on  image  at  t  =  1.20  =b  0.05  ns.  For  an  impact  pres¬ 
sure  ~  1  Gbar,  the  shock  speed  in  graphite  is  expected 
to  be  ~  250  fi m/ns  [12].  For  a  15  gm  thick  graphite,  this 
implies  shock  breakout  approximately  60  ps  after  the  col¬ 
lision  (assuming  a  steady  shock),  in  reasonable  agreement 
with  the  observation.  The  width  of  the  emission  break¬ 
out  is  approximately  200  gm  FWHM  at  the  time  of  the 
shock  breakout  and  300  gm  FWHM  at  the  time  of  the 
peak  emission  (110  ps  later),  indicating  that  the  width  of 
the  flyer  is  comparable  to  the  flat-top  of  the  laser  profile, 
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FIG.  6:  (Color  online)  Streaked  image  of  visible  light  emission 
from  the  rear  surface  of  the  second  foil.  Shock  breakout  occurs 
at  t  =  1.24  ns,  consistent  with  the  impact  time  from  the 
side-on  imaging.  The  spatial  extent  of  the  initial  breakout  is 
consistent  with  the  width  of  the  laser  hat-top. 


as  expected.  Significantly  the  shock  breakout  is  contin¬ 
uous  spatially,  indicating  that  the  target  integrity  was 
maintained  up  to  the  collision. 

Shots  with  significantly  thinner  foils  (<  10  fim)  allowed 
to  accelerate  over  longer  distances  («  800  gm)  show  clear 
signs  of  going  under  dense,  such  as  trajectory  breakup 
concurrent  with  the  appearance  of  self-emission  from  the 
laser  shining  through  to  the  second  foil. 


B.  Diagnosing  the  impact  via  neutron  production 

If  the  accelerated  foil  is  deuterated,  a  readily  mea¬ 
surable  neutron  yield  is  observed  on  impact.  Five  neu¬ 
tron  detectors  spaced  at  various  angles  around  the  target 
chamber  were  utilized.  Two  of  the  detectors  were  abso¬ 
lutely  calibrated  at  the  Gekko  XII  facility  in  Osaka  Uni¬ 
versity;  the  other  three  were  cross-calibrated  to  them.  A 
sample  neutron  signal  and  fit  are  shown  in  Fig.  7  for  a 
12  fim  CD  foil  initially  700  fim  from  the  graphite.  No 
significant  anisotropy  in  neutron  yield  was  found,  sup¬ 
porting  the  expectation  that  the  neutrons  are  thermonu¬ 
clear.  The  arrival  time  of  the  neutrons  is  as  expected 
for  2.45  MeV  DD  neutron  energy.  The  yield  on  this  shot 
was  measured  to  be  (1.3  ±  0.3)  x  106  neutrons.  For  an 
accelerated  CD  foil  when  the  second  foil  is  not  present, 
on  the  other  hand  (no  impact  takes  place),  the  yield  is 
measured  to  be  only  ~  5  x  104  neutrons  -  much  smaller, 
as  expected,  with  DD  reactions  taking  place  only  in  the 
laser  illuminated  low  density  corona. 

Ion  temperature  is  inferred  from  the  time-of-flight 
spread  in  the  neutron  signal  [13]  such  as  in  Fig.  7  as 
follows.  The  neutron  signal  is  fitted  by  a  convolution  of 
the  detector  response  with  the  Gaussian  spread  in  the 
neutron  arrival  times  due  to  the  ion  temperature.  The 
detector  response  is  a  convolution  of  an  exponential  de¬ 
cay  of  the  scintillator,  a  square  pulse  due  to  finite  thick¬ 
ness  of  the  scintillator,  and  the  impulse  response  of  the 
photomultiplier/cable/oscilloscope  system.  The  fitting 
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FIG.  7:  (Color  online)  Neutron  signal  (black)  and  convolu¬ 
tion  fit  (red)  for  a  12  /un  thick  CD  foil,  t  —  0,  as  before, 
corresponds  to  the  half-rise  of  the  laser  intensity  on  target; 
the  detector-target  distance  is  2.70  m;  given  that  the  impact 
occurred  at  t «  2  ns,  the  2.45  MeV  neutrons  are  expected  to 
arrive  at  t  =  127  ns,  consistent  with  the  time  of  the  rise  of 
the  signal. 


function  thus  has  the  parameters  of  the  Gaussian  (width, 
amplitude,  and  peak  location)  and  the  decay  time  of  the 
scintillator.  The  integral  of  the  Gaussian  is  proportional 
to  the  neutron  yield,  the  peak  location  gives  the  mean 
neutron  arrival  time,  and  the  FWHM  is  related  to  the 
ion  temperature:  Ti  oc  (FWHM)2  [14].  The  error  bar  on 
the  temperature  is  estimated  using  analysis  of  Ref.  15. 

An  interesting  result  was  seen  when  a  higher  density 
material  (gold)  was  used  in  place  of  graphite  second  foil 
in  some  of  the  shots.  The  neutron  yield  did  not  increase 
as  would  be  expected  at  first,  but  rather  was  an  order 
of  magnitude  lower  than  that  with  the  graphite  second 
foil.  A  likely  explanation  is  that  a  thin  surface  layer  of 
the  gold  is  pre-expanded  by  x-rays  from  the  laser  corona 
and  cushions  the  impact,  absorbing  a  significant  fraction 
of  the  flyer  kinetic  energy.  This  is  evident  in  the  side-on 
backlighting,  as  seen  in  Fig.  8.  Gold  from  the  second  foil 
is  clearly  expanding  toward  the  first  foil  during  the  tar¬ 
get  acceleration,  much  more  so  than  the  graphite  (Fig.  3). 
The  difference  is  likely  due  to  the  high  Z  of  Au,  which 
would  cause  it  to  absorb  the  x-rays  in  a  significantly  thin¬ 
ner  layer  than  graphite. 

IV.  DISCUSSION 
A.  Collision  hydrodynamics 

Useful  information  about  the  state  of  the  flyer  just  be¬ 
fore  impact  and  the  conditions  reached  upon  impact  can 
be  obtained  from  the  measured  neutron  yield  and  ion 
temperature.  In  order  to  interpret  these  results  analyti¬ 
cally,  the  hydrodynamics  of  the  collision  need  to  be  ex- 
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FIG.  8:  Significant  pre-expansion  is  observed  when  gold  is 
used  instead  of  graphite  as  the  second  foil.  Dashed  line  marks 
the  front  edge  of  material  expanding  from  the  second  foil. 


amined. 

For  a  CD  flyer  impacting  a  rigid  wall,  neglecting  ra¬ 
diation  and  conduction  losses,  the  kinetic  energy  of  the 
flyer  would  be  fully  converted  (less  ionization  energy)  into 
thermal,  heating  the  flyer  to  a  temperature  given  by  [16]: 

r  =  538  (urb:)^0077  keV-  (1) 

where  Vo  is  the  flyer  velocity  just  prior  to  the  impact. 

When  the  second  foil  is  not  a  rigid  wall  but  is  compress¬ 
ible,  some  of  the  kinetic  energy  of  the  flyer  is  imparted 
to  the  second  foil,  generating  a  moving  interface  between 
the  two  and  driving  a  shock  into  both.  This  is  illustrated 
in  Fig.  9,  which  shows  an  analytic  ID  hydrodynamics  x-t 
diagram  of  density  and  temperature  for  a  CD  foil  colli¬ 
sion  with  a  compressible  wall  (in  this  example  Cu).  In 
this  case,  Eq.  1  remains  valid  in  the  reference  frame  of 
the  contact  interface  of  the  foils.  The  temperature  ob¬ 
tained  in  the  collision  is  reduced  in  that  Vo  in  the  above  is 
substituted  by  an  effective  velocity  Vejj  —  Vo  —  Vinterface • 

Back-of-the-envelope  estimates  of  accelerated  foil  den¬ 
sity  and  unablated  mass  can  be  obtained  as  follows.  Tem¬ 
perature  T  and  neutron  yield  Yn  are  obtained  from  nTOF 
detectors  and  Vo  from  side-on  trajectory.  Veff  can  be  ob¬ 
tained  from  T  by  solving: 

r  = 538  (iftIvC  0077  keV>  (2) 

The  interface  velocity  Vinterface  is  then  obtained  from 
Vinterface  =  Vo  —  Veff  .  Equation  of  state  (EOS)  of  the 


FIG.  9:  Analytic  hydrodynamics  x-t  diagram  of  density  (left) 
and  temperature  (right)  for  a  CD  foil  collision  with  a  com¬ 
pressible  wall  (in  this  case  a  Cu  foil),  pen  —  0.8  g/cm3  and 
v  =  500  km/s.  Pressure  at  impact  is  1.55  Gbar,  shock  speed 
in  Cu  is  148  km/s,  Cu  compression  ratio  is  4.9,  and  the  inter¬ 
face  velocity  is  118  km/s,  reducing  the  effective  CD  velocity 
at  impact  to  382  km/s. 


second  foil  can  be  used  to  infer  the  collision  pressure  P 
knowing  the  interface  velocity.  The  flyer  density  prior  to 
the  collision  po  can  then  be  solved  for  from  the  expression 
for  P  obtained  assuming  ideal  gas  EOS  for  the  flyer: 

P=\p*Veff  (3) 

The  unablated  CD  areal  mass  can  be  computed  from 
the  neutron  reactivity  (ctv)dd  (depends  on  T  only)  and 
the  neutron  yield  Yn  by  solving  for  the  deuterium  areal 
number  density  Nd' 

V  3/  \  ANd2 

Yn  =  -^{ov)DD  y—  ,  (4) 

where  A  is  the  area  of  the  flyer. 

This  back-of-the-envelope  estimate  is  applied  to  mea¬ 
surements  from  a  particular  shot  in  the  table  shown  in 
Fig.  10.  Because  of  a  significant  error  bar  on  the  mea¬ 
sured  temperature  (2.7T0.4  keV),  the  inferred  quantities 
in  the  table  are  calculated  for  the  lower,  middle,  and  up¬ 
per  values  of  temperature  based  on  this  error  bar.  In 
addition,  for  a  particular  temperature,  the  inferred  un¬ 
ablated  mass  has  an  uncertainty  due  to  the  error  bar  on 
the  neutron  yield,  while  the  inferred  collision  pressure 
and  flyer  density  have  an  uncertainty  due  to  the  error 
bar  on  the  flyer  speed.  These  uncertainties  are  estimated 
by  error  propagation.  As  can  be  seen  in  the  table,  the 
temperature  uncertainty  has  the  dominant  effect  on  the 
inferred  values.  To  within  a  factor  of  2  however,  the 
inferred  collision  pressure  is  ^  1  Gbar,  inferred  flyer  den¬ 
sity  prior  to  impact  is  ~  l/10x  solid,  and  approximately 
1  pm  equivalent  of  solid  CD  remained  unablated. 

The  simple  analytic  estimates  above  neglect  a  lot  of 
physics  that  can  be  addressed  in  a  comprehensive  sim¬ 
ulation  of  the  collision,  such  as  realistic  density  profiles 
for  the  flyer  and  the  second  foil,  thermal  conduction  and 
radiation,  and  differences  in  ion  and  electron  tempera¬ 
tures.  For  example,  a  recent  analysis  of  such  planar  col¬ 
lisions  by  Gus’kov  et  al.  [17]  indicates  that  there  may  be 
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Measured  quantities 

Inferred  quantities 

T 

(keV) 

Neutron 

Yield 

(xIO6) 

Flyer 

speed 

(km/s) 

Effective 

speed 

(km/s) 

DD  neutron 
reactivity 
(cm3/s) 

Collision 

pressure 

(Gbar) 

Density  of 
flyer  (g/cc) 

Unablated  CD 
(/ym  solid  CD) 

2.3 

1.3  ±0.3 

910  ±30 

665 

4.8E-21 

1.54  ±0.38 

0.26  ±0.06 

1.39  ±0.16 

2.7 

1.3  ±0.3 

910  ±30 

718 

9.4E-21 

0.94  ±  0.30 

0.14  ±0.04 

1.03  ±0.12 

3.1 

1.3  ±0.3 

910  ±  30 

768 

1 .6E-20 

0.52  ±  0.22 

0.07  ±0.03 

0.81  ±0.09 

FIG.  10:  “Back-of-the-envelope”  estimates  of  flyer  density 
and  remaining  mass  based  on  measured  temperature,  neutron 
yield  and  target  velocity. 


insufficient  time  for  electron  and  ion  temperature  equili¬ 
bration.  If  this  were  the  case  it  would  imply  that  a  lower 
Veff  corresponds  to  the  measured  ion  temperature.  The 
inferred  unablated  mass  is  only  weakly  affected  by  this 
however,  which  can  be  seen  as  follows.  Generalizing  Eq.  2 
for  Te  ^  Ti  and  neglecting  the  ionization  loss  (0.077  keV) 
we  get: 

3  {kTi  +  ZkTe ) 

Mmn  ’  ^  ] 

where  Z  =  3.5  is  the  mean  charge  and  M  —  7  the  mean 
atomic  weight  of  CD  and  mn  is  the  atomic  mass  unit. 
From  Eq.  4,  for  a  given  ion  temperature  and  neutron 
yield,  N&  oc  \fVe^  oc  (Ti  +  ZTe)i.  Thus,  for  exam¬ 
ple,  if  Te  =  1/2  Ti,  Nd  would  only  be  reduced  by  12%. 
Lower  Veff  would  likely  imply  however,  that  the  second 
foil  is  pre-expanded,  cushioning  the  impact.  Future  ex¬ 
periments  will  aim  to  measure  Te  spectroscopically  using 
a  high-Z  dopant  in  the  flyer. 


B.  Simulations 

A  significant  benefit  of  the  experiments  like  those  dis¬ 
cussed  here  is  that  they  provide  a  useful  platform  for 
exercising  and  benchmarking  radiation  hydrodynamics 
codes.  Some  initial  results  of  simulations  of  the  flyer 
acceleration  are  presented  in  this  section.  Comprehen¬ 
sive  simulations  of  the  collision  will  be  the  subject  of  a 
separate  article. 

An  important  experimental  result  is  that  the  target 
acceleration  was  sufficiently  robust  against  the  Raleigh- 
Taylor  instability  to  maintain  its  integrity  over  a  dis¬ 
tance  ~  60  x  its  thickness.  This  result  can  be  examined 
first  by  applying  analytic  RTI  linear  growth  formula  to  a 
simple  ID  simulation  of  the  acceleration.  Fig.  11  shows 
the  time-integrated  amplification  of  linear  perturbations 
computed  by  post-processing  a  ID  simulation  of  a  12  gm 
CD  foil  driven  by  the  pulse  in  Fig.  2.  The  amplifica¬ 
tion  is  computed  by  integrating  the  growth  rate  given  by 
the  modified  Takabe  formula  [18].  As  can  be  seen  in  the 
figure,  the  maximum  amplification  factor  is  only  about 
lOOx.  Given  the  ISI  laser  smoothing  and  a  smooth  foil 
(foils  used  in  the  experiment  had  0.4  nm  rms  roughness), 
this  implies  perturbation  growth  to  significantly  less  than 
1  f± m  amplitude  -  insufficient  to  break  up  the  foil. 


FIG.  Ill  Time-integrated  amplification  of  linear  perturba¬ 
tions  (post-processing  of  ID  simulation).  The  parameters  of 
the  modified  Takabe  formula  shown  are  taken  from  the  simu¬ 
lation. 


In  order  to  examine  stability  as  well  as  the  effects  of  a 
finite  spot  size,  a  high  resolution  2D  axisymmetric  sim¬ 
ulation  using  NRL’s  FASTRAD3D  radiation  hydrocode 
has  been  performed.  Fig.  12  shows  a  plot  of  mass  density 
from  the  simulation  at  a  point  in  time  when  the  target  has 
traveled  just  over  550  /xm  and  reached  1000  km/s.  The 
target  is  a  12  /x m  CD  foil  driven  by  the  pulseshape  and 
spatial  profile  used  in  the  experiment  (Fig.  2);  the  simu¬ 
lation  mesh  size  is  0.1  /xm  along  the  laser  axis  and  3.5  /x m 
perpendicular  to  it.  Although  the  finite-spot-size  effect  is 
prominent,  the  target  has  maintained  its  integrity.  The 
peak  density  is  0.12  g/cm3,  not  far  from  the  estimates  ob¬ 
tained  in  Sec.  IV  A  above.  The  “fiat-top”  portion  of  the 
target  is  200-300  /xm,  consistent  with  the  shock  breakout 
measurements  (Fig.  6). 

In  order  to  examine  the  target  trajectory,  on- axis  val¬ 
ues  of  peak  density,  velocity,  and  distance  traveled  are 
taken  from  the  2D  simulation  above  and  are  shown  in 
Fig.  13.  As  can  be  seen  in  the  figure,  peak  density  de¬ 
cays  after  the  initial  shocks,  but  stabilizes  during  the 
later  portion  of  the  trajectory.  The  acceleration  increases 
toward  the  end,  similar  to  what  is  seen  in  the  experiment 
(Fig.  5). 

For  a  more  direct  comparison  with  the  shot  shown  in 
Fig.  4,  figure  14  shows  the  on- axis  lineouts  of  a  2D  ax¬ 
isymmetric  simulation  of  a  10.5  /xm  CH  target.  In  addi¬ 
tion  to  the  trajectory  and  velocity  of  the  peak  of  mass 
density,  the  plot  shows  the  trajectory  and  the  velocity  of 
the  edge  of  the  x-ray  emission  in  the  1.25-2.5  keV  band. 
This  x-ray  energy  band  is  centered  on  the  1.86  keV  x-ray 
energy  imaged  in  the  experiment  (Fig.  4).  Both  trajecto¬ 
ries  overlap  to  within  the  line  thickness,  while  the  veloc¬ 
ities  differ  significantly  only  above  1100  km/s,  where  the 
density  is  dropping  precipitously.  Increase  in  accelera¬ 
tion  prior  to  that  is  similar  to  that  seen  in  the  experiment 
(Fig.  5). 
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1.41  ns;  V=980  km/s;  p  =0.12  g/cm3 

rmax  3 


-700  -600  -500  -400  -300 

z  (urn) 


FIG.  12:  (Color  online)  Mass  density  from  an  axisymmetric 
2D  simulation  of  a  12  /rm  CD  foil  after  the  target  has  traveled 
about  550  /mi.  The  z-axis  at  the  bottom  of  the  figure  is  the 
axis  of  symmetry  with  the  laser  shining  from  the  right;  the 
simulation  mesh  size  is  0.1  /mi  along  the  laser  axis  and  3.5  /mi 
perpendicular  to  it.  Peak  density  is  0.12  g/cm3,  velocity  is 
980  km/s,  and  “flat-top”  portion  of  the  target  is  200-300  /mi. 


portunities.  Utilizing  a  KrF  laser  allowed  us  to  achieve 
this  with  as  little  as  1  kJ  of  energy.  Higher  velocities 
should  be  possible:  higher  laser  energy  would  allow  use 
of  thicker  targets  that  should  reach  higher  velocities  still. 

Impact  of  the  highly  accelerated  targets  results  in 
~  Gbar  pressures  and  heating  to  thermonuclear  tempera¬ 
tures.  Neutron  production  gives  an  effective  diagnostic  of 
the  state  upon  impact.  However,  higher  densities  of  the 
accelerated  foils  are  needed  for  impact  ignition.  Future 
experiments  will  attempt  to  obtain  these  velocities  with 
less  preheat,  which  would  allow  higher  in-flight  densities. 
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FIG.  13:  (Color  online)  Peak  density,  target  velocity,  and 
distance  traveled  on  the  axis  of  the  2D  simulation  of  Fig.  12 
of  a  12  /mi  CD  foil. 


V.  CONCLUSIONS 

Target  velocities  on  the  order  of  1000  km/s  have  been 
routinely  achieved  for  the  first  time,  opening  up  new  op- 


FIG.  14:  (Color  online)  Peak  density,  target  velocity,  and 
distance  traveled  on  the  axis  of  a  2D  axisymmetric  simulation 
of  a  10.5  /mi  CH  foil.  Increase  in  acceleration  at  the  end  is 
similar  to  that  seen  in  the  experiment  (Fig.  5).  In  addition 
to  the  trajectory  and  velocity  of  the  density  maximum,  the 
trajectory  and  velocity  of  the  edge  of  the  x-ray  emission  in 
the  1.25-2.5  keV  band  are  plotted.  Both  trajectories  overlap 
to  within  the  line  thickness  (labelled  “distance”),  while  the 
velocities  differ  significantly  only  above  1100  km/s. 
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